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ABSTRACT 

We perform a spectral analysis of a sample of 11 medium redshift (1.5 < z < 2.2) 
quasars. Our sample all have optical spectra from the SDSS, infrared spectra from GNIRS and 
TSPEC, and X-ray spectra from XMM-Newton. We first analyse the Balmer broad emission 
line profiles which are shifted into the IR spectra to constrain black hole masses. Then we 
fit an energy-conserving, three component accretion model of the broadband spectral energy 
distribution (SED) to our multiwavelength data. Five out of the 11 quasars show evidence of 
an SED peak, allowing us to constrain their bolometric luminosity from these models and 
estimate their mass accretion rates. Based on our limited sample, we suggest that estimating 
bolometric luminosities from T- |(|0A and L 2-10 keV may be unreliable, as has been also noted 
for a low-redshift. X-ray selected AGN sample. 

Key words: black hole physics; accretion discs; quasars: supermassive black holes; galaxies: 
active, high-redshift 


1 INTRODUCTION 
1.1 Background 

Multi-wavelength studies of quasars are an important means to fur¬ 
ther our understanding of these objects. Interpreting the quasar 
spectral energy distribution (SED) remains a challenge in astro¬ 
physics, although the last few decades have seen significant ad¬ 
vances in the field (e.g. |Ward et ah||1987| |Elvis et ah||1994| |Va-| 
|sudevan & Fabian 2007}. It is now widely accepted that accretion 

of gas onto a central supermassive black hole (SMBH) is the ulti¬ 
mate source powering these extremely luminous objects. The SED 
contains clues to the geometry and properties of the matter in the 
regions close to the black hole (BH). Better understanding of the 
SEDs also has consequences for cosmology, as there is strong evi¬ 
dence that galaxy formation is influenced by the quasar (hereafter 
referred to as an active galactic nucleus, AGN) through the process 
of feedback (e.g. |McCarthy et al.|2010| ). 

The simplest interpretation of the unified model of AGN pro¬ 
poses that most differences in the observed SED properties can 
be attributed to differences in the orientation of the accretion flow 
and line-of-sight absorption. Those AGN characterised by a direct 
sightline towards the central engine are classed as Type I, and those 
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with optically thick material along the line-of-sight are classified 
as Type II l )Antonuccill993| l. In practice, however, this simple pic¬ 
ture is incomplete, and other physical differences in the BH itself 
(e.g. spin) and the properties of the infalling matter influence the 
observed SED ( [Boroson & Green|1992[|Done et al. [2012k 

This situation has motivated us to produce new broadband 
SED models, based on theoretical considerations and on empiri¬ 
cal studies of accreting stellar-mass BHs found in X-ray binaries. 
Combining AGN SED models with representations of galactic and 
extragalactic extinction, both via dust and photoelectic absorption, 
enables fitting of the multi-wavelength data from an AGN, and re¬ 
covery of the intrinsic SED, which relates directly to the properties 
of the BH and the material it accretes. Previous studies by Jin et 
al. (2012a, b, c), hereafter J12a, b, c (collectively J12) have suc¬ 
cessfully employed this technique. The peak of the energy output 
occurs at far-ultraviolet/ultra soft X-ray energies in the majority 
of these low redshift AGN. These energy ranges are mostly un¬ 
observable due to photoelectric absorption by neutral Hydrogen in 
the inter-galactic medium along the line-of-sight to the source (the 
so-called Gunn-Peterson trough, [Gunn & Peterson|1965| >, and that 
intrinsic to the Milky Way. 

In order to overcome this restriction we consider higher red¬ 
shift AGN, in which it is expected that the peak of the SED would 
be shifted into the observable optical/near-UV energy range. This 
occurs for two reasons; the first is simply the redshift, and the sec- 
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Figure 1. An SED characteristic of a typical disc-dominated AGN. Colours 
correspond to the regions shown in the AGN schematic (top), with red rep¬ 
resenting the AD, green the SX and blue the PLT from the corona. The 
grey dotted line shows the resultant observed SED when it has been attenu¬ 
ated by typical galactic and extra-galactic extinction and absorption. Figure 
adapted from Done et al. (2012). 


ond is because these more distant, luminous AGN contain more 
massive BHs which have cooler accretion discs peaking at lower 
energies fMcLure & Dunlop|2004[|Done et al.|2012| l. 

In this first paper in a series, we present the sample and dis¬ 
cuss the selection process and the data that has been assembled. We 
then present model fits from rest-frame optical to hard X-rays, and 
explore the modifying effects of extinction and the presence of a 
stellar component in the host galaxy. In the next paper (Collinson 
et al. (in prep), hereafter Paper II), we will investigate the parameter 
space further, including the toroidal dust component. In Paper III, 
we will apply these findings to the analysis of a larger, statistically 
significant sample. 


1.2 A Refined AGN SED Model 

We use a multi-component model based on studies of black hole 
binaries (BHBs) and nearby Narrow-line Seyfert Is (NLSls). This 
model is described in |Done et al.| < ]2012) . It is characterised by 
three principal components: an accretion disc (AD), a power-law 
tail (PLT) and a soft X-ray excess (SX). A schematic SED diagram 
is shown in Fig. |T] The AD is modelled as a relativistic, geomet¬ 
rically thin, optically thick disc, with each radius in the disc radi¬ 
ating as a blackbody (e.g. |Shakura & Sunyaev||1973) . Our latest 
disc models include a colour correction (/ co i) to account for the 
fact that the disc is not completely thermalised at all radii. There 
is also a large contribution from the PLT at high energies, arising 
from inverse Compton scattering of AD photons by a hot, optically 
thin corona (Zdziarski et al. 1995}. 

The origin of the SX is less well understood, and is the sub¬ 
ject of debate. Some postulate that the SX is produced by reflec¬ 
tion of hard coronal X-rays off the AD (e.g. Crummy et al.|2006[ 


|Fabian et al.||2009| |Zoghbi et al.||201l) , whilst others attribute it 
to the presence of intervening matter, which complicates the ob¬ 
served emission (e.g. jMiller et al.|2010j >, and that a warm, optically 
thick Comptonised component of the inner AD better describes the 
observations (e.g. |Alston et al.|2014||Gardner & Done|2014}. |jin| 
|et al. | (2013) and |Matt et al. 1(2014) used long observations with 
XMM-Newton and NuSTAR to test models of the SX in two dif¬ 
ferent AGN, and also found Comptonisation of inner disc photons 
to better describe their observations. Whilst the model we adopt 
in this paper assumes the latter origin for the SX, in practice lim¬ 
itations in the quality and energy coverage of our data mean that 
uncertainties in the origin of the SX cannot be unambiguously re¬ 
solved. Additional factors arising from absorption and reflection 
could be included. The interplay between the model components is 
complex, but by fitting such physically motivated models to data, 
we can infer information about the accretion flow properties (e.g. 
|Elvis et al,| 1994] |E1 vis et al.|2012) . 

Importantly, the Done et al. | (2Q12) AGN SED model we use 
- OPTXAGNF - applies the constraint of energy conservation, as 
energy output is dependent on the amount of matter accreting on 
to the BH. This model has so far been tested on a moderate num¬ 
ber of AGN (e.g. J12, |Done et al.||20T3l|Matt et al.|[20l4} . It has 
been designed for implementation into NASA’s High-Energy As¬ 
trophysics Science Archive Research Center (HEASARC) XSPEC 
spectral fitting package. 


1.3 Previous Work 

Initial studies of AGN focused primarily on ultraviolet (UV), op¬ 
tical and infrared (IR) spectra (e.g. |Wills et al,|1985[|Zheng et al.| 
|1997| >, but with improvements in ground-based instrumentation and 
new satellites, it is now routinely possible to study AGN samples 
across multiple wavelength bands. 

|Puchnarewicz et aT7| ( | 1992| > presented a study of optical and X- 
ray data of 53 AGN with ultra-soft X-ray excesses, and found a bias 
in their emission line profiles towards narrow linewidths, which 
had implications for the position and size of the broad line region 
(BLR). |Grupe et al.](1998) and |Grupe et al. (1999) con firmed this 
result and also reinforced the findings of |Walter & PTTTk]ff993> that 
there is a Big Blue Bump (BBB) from optical to X-ray spectra. 

More recently, J12 presented a medium-sized SED modelling 
study of 51 AGN. In their study, they assembled optical and X-ray 
spectra from the Sloan Digital Sky Survey (SDSS) and the Euro¬ 
pean Space Agency (ESA) XMM-Newton, respectively. This was 
supplemented, when available, by photometric UV data from the 
XMM Optical Monitor (OM). These data ranges are not contiguous 
and leave the SED peak, unobservable due to the aforementioned 
absorption, devoid of data. Whilst SDSS and XMM OM data con¬ 
strain only the optical edge of the AD, XMM EPIC data only lies in 
the energy range of the SX and PLT. Thus the AD peak (and hence 
SED peak in the disc-dominated objects) was unsampled and in¬ 
formation as to its position and shape had to be inferred from the 
fitted models. Since the AD peak originates from the innermost part 
of the accretion flow, it contains key information about the BH spin 
and mass accretion rate, m, and hence the bolometric luminosity of 
the AGN. 

BHs have ‘no-hair’, meaning they are characterised by just 
three properties - mass, spin and charge. Much study has been 
devoted to measuring the mass of a BH (see §1.4) , and charge is 
negligible in an astrophysical context. There is currently much in¬ 
terest in constraining spin, which can be estimated if some mea¬ 
sure of the radius of innermost stable circular orbit (ri sco ) can be 
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made. In practice this is difficult; Fabian et al. (2009) and Risaliti 
|et al.| j2013) claim to make such measures from observations of the 
broad iron K emission line, but this is still controversial (e.g. |Millerj 
|& Turner|2013| >. 

|Davis & Laor| ( |20fT] l, |Done et al,| ( [20T3} , |Matt et al. | (20~i~4| i and 
|Capellupo et al.|f201 have all investigated the potential for using 
SED modelling to constrain BH spin. |Done et al.|( 2013 ) utilised the 
same energy conserving model discussed in j |1.2| to rule out a high- 
spin solution for the SMBH at the centre of PG1244+026. |Matt| 
|et al.||2014| l also used this model in their study of Ark 120, infer¬ 
ring an intermediate BH spin from simultaneous XMM-Newton and 
NuSTAR observations. Capellupo e t al.| |20151 fitted a thin accretion 
disc model ( jSlone & Netzer|2012 1 to the IR/optical continua of 30 
quasars in a similar redshift range to our sample. Using this model 
they inferred a range of BH spin parameters in their sample. 

Whilst previous works have successfully found evidence for 
the AD turnover by recovering the continuum from high-resolution 
UV spectra of high-redshift AGN, (e.g. |Zheng et al.||1997| |Finn| 
|et al.|2014] > no sizeable sample exists with X-ray spectra. Adding 
to this is the need for reliable BH mass estimates. 


1.4 Black Hole Mass Estimates 


Reverberation Mapping (RM), proposed by |Blandford_^^lcKee 


1982) and employed by e.g. Peterson et al. 1 2004), Bentz et al. 


(2009 and |Denney et al.| < |2010|) , amongst others, is a technique 


for accurately determining the mass of BHs in AGN. RM involves 
measuring the delay in time for variations in the central source in¬ 
tensity (continuum) to reach the broad line region (broad emission 
lines), and using this as a proxy for the light travel time across 
the BLR, which is well-correlated with the BH mass. RM requires 
extensive long-term monitoring programs, and is thus very expen¬ 
sive observationally, however, it has allowed the calibration of BH 
mass estimates from single-epoch observations of the broad emis¬ 
sion lines in AGN optical spectra (e.g. |Woo & Urry|2002, Greene 
& Ho 2005 Mat suoka et al.|2013fr . The abundance of suitable sin¬ 
gle epoch spectra means that this method has now been applied to a 
great many AGN, generally utilising the well-studied Balmer-series 
Hydrogen lines Ha and H f3. 

The problem with this technique arises in high-redshift AGN, 
where the Balmer lines are redshifted to infrared wavelengths, and 
rest-frame UV lines are shifted to the optical regime. Much effort 
has been devoted to the use of the rest-UV lines for the purpose of 
BH mass estimation, but the subject remains contentious. Vester- 
Igaard & Peterson (2006 ) studied the scaling relationships between 
mass estimates from H/3, Ha and the rest-UV lines CIV and Mg II, 
finding that CIV and Mg II offer viable alternatives to the Balmer 
lines for this purpose in a sample of low-redshift AGN. However, 
two studies of high-redshift AGN, |Netzer et al.|{2007fr and |Shen &] 
|Liu ] p0l2) , both disagree with this result, contending that the line 
profile of CIV is not suitable for mass estimates, but agreeing that 
Mg II shows reasonable correlation. Possible explanations put for¬ 
ward for this discrepancy in CIV are the presence of outflows that 
influence line profiles in high-redshift, high-Eddington-ratio AGN. 
It has also been claimed l |Denney et al.|201~3] > that high signal-to- 
noise ratio (S/N) spectra are required to derive reliable mass esti¬ 
mates from CIV. In any case, H/I remains the mass estimator of 
choice, i jWoo & Urry||20 02|> s imply by its virtue as the line best 
calibrated by RM. Greene & Ho| (2005| l propose an alternative in 
Ha which shows excellent correlation with H/?, but which, due to 
greater line strength, generally offers better S/N for the line analy- 


2 SAMPLE AND DATA ASSEMBLY 

2.1 Sample Selection 

In order to define an AGN sample in which the SED peak would be 
observable, we required objects with redshifts around z ~ 2. For 
SED modelling, we require data lying on both sides of the UV/soft 
X-ray absorption trough, as was the case in the J12 sample. We 
thus required optical spectra, available from the SDSS (~ 35% sky 
coverage) and X-ray data, available from XMM-Newton s 3XMM 
DR4 data catalogue (~ 2% sky coverage). In order to constrain 
physical parameters for the SED model we also required IR spectral 
data, so that BH mass estimates from the Balmer lines could be 
made. 

We started by searching the |Schneider et al.| i |201()f SDSS DR7 
QSO catalog for all AGN meeting the following criteria: 

• 1.49 < z < 1.61: This was so that Ha and H/I would lie in the 
NIR H and J bands respectively, and Mg II and CIV lines would 
be visible in the SDSS spectra, for comparison purposes. 

• Admass < 16.5: Since we required medium resolution IR 
spectra, we needed the objects to be suitably bright in the IR bands. 
Where K was unavailable in 2MASS, similar constraints were ap¬ 
plied to the J and H bands. There were 1797 matching AGN after 
this step. 

• XMM-Newton data: We also required X-ray spectra from 
ESA’s XMM-Newton instrument. 63 objects had matching observa¬ 
tions with 18 of these being bright enough to have had a spectrum 
extracted by the XMM pipeline (typically requiring > 200 counts). 

We were awarded observing time in Cycle 2013B to use 
the Gemini GNIRS instrument, sited at Mauna Kea Observatory, 
Hawai’i, to obtain high quality IR spectra in the J, H and K bands. 
In practice 9 of our objects were visible to GNIRS in 2013B, and 
we proposed to observe the 6 with the highest SDSS S/N. 

Simultaneously, we identified 4 objects in the |Shen & Liu| 
|20l2| sample with XMM X-ray data of varying quality, and IR 
spectral data from the ARC TripleSpec (TSPEC) instrument. One 
additional object with an XMM observation was selected from a 
publicly available, archival GNIRS 2004B dataset (project GS- 
2004B-Q16, PI Todd Boroson) that had not previously been re¬ 
duced. This resulted in our final sample of 11 objects, listed in Table 
□ 

It is plausible that our selection criteria introduces some bias 
into the sample, for instance, by choosing objects with a higher 
than average IR/optical or X-ray/optical ratio. However, our mean 
value of aox for Model 3 (see §4.3| ) is 1.61 ± 0.09, which is typi¬ 
cal of large X-ray selected samples (e.g.|Vignali et al.|2003||Lusso| 
|et al.|2010[ >. There may be an inherent bias in our sample owing 
to the flux threshold of the parent SDSS spectroscopic database. It 
is of course true that the sample must have X-ray data, which will 
naturally bias against some sub-classes of AGN, e.g. BAL quasars. 

2.2 Optical/IR Data Preparation 

The IR spectra for the four objects selected from the [Shen & Liu| 
( |2012| ) sample were kindly provided by Yue Shen (private commu¬ 
nication). 

The data resulting from Gemini cycle 2013B were reduced 
according to the guidelines provided on the Gemini website, us¬ 
ing IRAF V2.14, and the Gemini IRAF package Vl.12. The archival 
GNIRS 2004B object (J0118—0052) was reduced before the re¬ 
lease of Vl.12 and we thus used Vl.l 1.1 of the Gemini IRAF pack¬ 
age to reduce this object. 
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Table 1. The names, positions, SDSS pipeline redshifts and UV/IR data sources for the sample of 11 objects. The UV photometry is of limited use for our 
purposes (see j|4.2|. 


ID 

Name 

R.A. (J2000) 

Decl. (J2000) 

Common Name 

ZSDSS 

IR Source 

UV Phot a 

1 

J0041-0947 

00 41 49.64 

-09 47 05.0 

— 

1.629 

TSPEC 6 

N 

2 

J0043+0114 

00 43 15.08 

+01 14 45.6 

- 

1.563 

GNIRS 13B 

N 

3 

J0118-0052 

01 18 27.98 

-00 52 39.8 

QSO B0115-0108 

2.188 

GNIRS 04B 

- 

4 

J0157—0048 

01 57 33.87 

-00 48 24.4 

QSO JO 157-0048 

1.551 

TSPEC 6 

NM2W1 

5 

J0839+5754 

08 39 06.53 

+57 54 17.0 

3C 205 

1.534 

GNIRS 13B 

NW1 

6 

J1021+1315 

10 21 17.74 

+13 15 45.9 

- 

1.565 

GNIRS 13B 

W1 B 

7 

J1044+2128 

10 44 01.13 

+21 28 03.9 

- 

1.494 

GNIRS 13B 

N 

8 

J1240+4740 

12 40 06.70 

+47 40 03.3 

- 

1.561 

TSPEC 6 

NFU 

9 

J1350+2652 

13 50 23.68 

+26 52 43.1 

QSO B1348+2707 

1.624 

TSPEC 6 

NF 

10 

J2328+1500 

23 28 10.56 

+15 00 12.8 

- 

1.536 

GNIRS 13B 

N 

11 

J2332+0000 

23 32 28.21 

+00 00 32.8 

- 

1.604 

GNIRS 13B 

N 


°N = OT/Vgalex, F = TWgalex- M2 = UVM2 x mm qm.WI =U VWIxmm om, B =Bxmm om 
6 IR spectral data from TSPEC courtesy of Yue Shen 


Shen & Liu 


2012 . 


Table 2. Optical/IR spectral observation dates for our sample of 11 objects. S/N values are approximate. 


ID 

Name 

Survey 

MJD-Plate-Fibre 

Optical 
Obs. UT 

Avg. 

S/N 

IR 

Source 

IR 

Obs. UT 

Avg. 

S/N 

1 

J0041—0947 

SDSS 

52162-655-172 

2001-09-10 

36 

TSPEC 

2010-01-02 & 11-28 

7 

2 

J0043+0114 

SDSS 

51794-393-419 

2000-09-07 

18 

GNIRS 13B 

2013-08-16 

12 



BOSS 

55186-3589-707 

2009-12-21 

30 






BOSS 

55444-4222-902 

2010-09-05 

33 




3 

J0118-0052 

SDSS 

51789-398-211 

2000-09-02 

16 

GNIRS 04B 

2004-11-29 

9 

4 

J0157—0048 

SDSS 

51871-403-213 

2000-11-23 

17 

TSPEC 

2009-11-07 & 11-28 

3 



SDSS 

52179-701-294 

2001-09-27 

14 






BOSS 

55449-4233-152 

2010-09-10 

24 




5 

J0839+5754 

SDSS 

54425-1784-495 

2007-11-21 

32 

GNIRS 13B 

2013-10-27 

19 

6 

J1021 + 1315 

SDSS 

53062-1746-491 

2004-02-27 

14 

GNIRS 13B 

2014-03-21 

15 

7 

J1044+2128 

SDSS 

54097-2478-411 

2006-12-28 

18 

GNIRS 13B 

2014-03-20 

14 



BOSS 

56039-5874-970 

2012-04-22 

28 




8 

J1240+4740 

SDSS 

53089-1455-424 

2004-03-25 

16 

TSPEC 

2011-02-22 

8 

9 

J1350+2652 

SDSS 

53848-2114-105 

2006-04-23 

27 

TSPEC 

2011-02-22 

8 



BOSS 

56105-6006-260 

2012-06-27 

47 




10 

J2328+1500 

SDSS 

52238-746-463 

2001-11-25 

10 

GNIRS 13B 

2013-08-18 

11 

11 

J2332+0000 

SDSS 

51821-384-438 

2000-10-04 

13 

GNIRS 13B 

2013-08-19 

13 



SDSS 

52199-681-543 

2001-10-17 

13 






SDSS 

52525-682-355 

2002-09-08 

11 





All Gemini spectra were created by GNIRS in cross-dispersed 
(“XD") mode, and we were able to recover orders 3-8 for the 2013B 
objects. We recovered orders 3-6 in the 2004B object, due to a more 
limited range of flat fields provided. All stages of the reduction 
were visually inspected, to ensure no errors had occurred. 

Telluric stellar spectra were provided by Gemini. The purpose 
of these spectra are to provide a well-defined reference spectrum 
that may be used to correct object spectra for non-constant sky 
transmission across the infrared wavelength range. To correct for 
telluric features, we were not able to use the XTELLCOR routine 
discussed in |Vacca et al.|{2003| l as not all of the telluric stars ob¬ 
served were of spectral class AOV. Hence we corrected for telluric 
features in the following way. 


Hydrogen absorption features in the telluric spectrum were 
first removed in IRAF, as these features are intrinsic to the star it¬ 
self. Lorentzian absorption line profiles were assumed, as these fit 
the features better than Gaussian or Voigt profiles. The assumed 
template for the telluric star was a blackbody of a characteristic 
temperature dependent on the telluric stellar class, and flux nor¬ 
malised to the 2MASS J magnitude. The extracted object spectrum 
was than multiplied by the ratio of the assumed telluric template to 
the extracted telluric spectrum, producing a very satisfactory cor¬ 
rection to the variable atmospheric transmission. 

Two factors may affect the relative normalisations of the IR 
and optical spectra, observed as they are at different epochs. Vari¬ 
ability of the AGN is one possible origin for such a change, and 


© 2014 RAS, MNRAS 000 .ftp!] 























Reaching the Peak of the Quasar Spectral Energy Distribution 1: Observations and Models 5 


will be discussed more in j (5.3| The other factor is the accuracy 
of the flux calibration; if the optical or IR spectra are flux cal¬ 
ibrated poorly (e.g. due to seeing problems or aperture effects), 
the resulting normalisation of the spectra will be incorrect. For the 
SDSS/BOSS spectra, we will assume that the flux calibration is 
precise. 

The objects we reduced from Gemini were flux calibrated in 
the telluric correction step discussed above, which makes the as¬ 
sumption that the telluric star is not itself variable, and that weather 
conditions do not change between the object and the telluric obser¬ 
vations. The Shen TSPEC objects, on the other hand, are all flux- 
normalised to the 2MASS H band magnitudes. 

All of the objects in our sample have either UKIDSS or 
2MASS (or both) photometry in each of J, H, and K bands. We 
first compare the flux of each photometry point with the average 
flux in the IR spectrum across the effective photometric bandpass, 
to estimate the percentage difference in flux density between these 
two values. As an additional check, we can also compare how well 
the blue wing of each IR spectrum fits a power-law extrapolated 
redwards of the SDSS spectrum, or, if available, the region of spec¬ 
tral overlap between optical/IR spectra. This highlights offsets of 
the IR spectrum relative to the optical spectrum. In doing this, we 
find that all but three of our objects show evidence for less than 
10% difference between the IR spectrum and optical spectrum/IR 
photometry. 

The three objects showing greater than 10% variation be¬ 
tween optical and IR spectra are 10041—0947, J1044+2128 and 
J2328+1500. 

In J1044+2128 and J2328+1500, we normalise to the A'-band 
photometry, as there is evidence for some cloud cover at the time of 
observation that could affect the flux calibration. In J0041—0947, 
the IR spectrum was normalised to 2MASS H by Shen & Liu. How¬ 
ever, it lies below the level predicted from the SDSS spectrum, 
possibly due to variability. We thus make our mass estimate from 
the spectrum as flux-normalised by Shen & Liu, but normalise the 
spectrum to the optical level for the SED fit, as whether the differ¬ 
ence is due to flux calibration error in either spectrum, or variability, 
we require agreement to fit the SED shape. 

Finally, in all objects we corrected both the IR and optical 
spectra for extinction by the Milky Way, using the dust maps of 
Schlegel et al.| (|1998| and the extinction law of |Cardelli et al.| 
l |1989| i. 

All of our optical data come from the SDSS. Five of the ob¬ 
jects had been observed multiple times in the SDSS spectroscopic 
survey, or reobserved by BOSS, the follow-up to SDSS-III. The 
BOSS survey offers greater wavelength coverage and improved 
S/N for the observed objects, and hence the BOSS spectra are our 
preferred source of optical data. In the absence of BOSS spectra, 
we use the SDSS spectrum that best aligns with the IR spectrum. 

The optical and IR spectral observation dates are tabulated in 
Table[2] We provide estimates of the spectral S/N ratio by estimat¬ 
ing and averaging the S/N in 10000 random, 50-pixel subsamples 
from each spectrum. The optical spectra adopted for the SED fit¬ 
ting are highlighted in bold. A greater discussion of variability in 
general is given in §5.3| 

2.3 UY and X-ray Data 

We extracted all X-ray data from the XMM Science Archive, 
through the HEASARC. For each object, we obtained spectra 
for all available observations from the EPIC MOS1, MOS2 and 
PN cameras, to maximise the number of counts. Each spectrum 


comes in three parts: a source spectrum, a background spectrum, 
and an ancillary response function. These are supplemented by an 
instrument-specific canned response matrix, downloaded from the 
XMM-Newton EPIC Response Files Page. These are all provided in 
a file format readable by XSPEC. 

For two of the objects in our sample for which IR spectra were 
pre-existing, the X-ray spectrum had not been extracted from the 
EPIC source image, due to lack of counts. However, the XMM- 
Newton Serendipitous Source Catalog (3XMM DR4) lists, for all 
detections, the X-ray flux in five bands: 0.2-0.5 keV, 0.5-1.0 keV, 
1.0-2.0 keV, 2.0-4.5 keV, 4.5-12.0 keV, and we thus used these 
points in the SED fitting. 

XMM-Newton also has an optical monitor (OM), which can 
obtain UV photometry simultaneously to the EPIC observation. 
Targets that are too far off-axis are not covered, so many of our 
sources that were observed serendipitously by EPIC do not have 
OM data. 

Finally, we also searched the Galaxy Evolution Explorer 
( GALEX) database for UV photometry. All sources have some data 
in the GALEX All-sky Imaging Survey (AIS). 


3 BLACK HOLE MASSES 

The first step in analysing data was to make BH mass estimates 
from the IR spectra. Discussion of a full spectral decomposition 
will be made in Paper II, but relies on the results of our broadband 
SED fits, so at this stage we will use localised decomposition of the 
B aimer line region, which is redshifted into our IR spectra. 

3.1 Optical/IR Spectral Fitting 

We first perform a continuum subtraction of the region under the 
Balmer lines. Whilst the underlying, featureless continuum is in 
theory best described by an accretion disc, on localised scales it 
may approximate a power-law continuum. As discussed in e.g. 
|Vanden Berk et ah] ( (2001 [ ) the continuum under the Balmer lines 
is more accurately represented by a broken power law, so in our 
Balmer region continuum fit we also employ a broken power-law. 
We use, for the power-law, the standard form: 

F(X) = Ci (A/5100 A) -C2 (1) 

where Ci is the normalisation and C 2 is the power-law slope, and 
in our case, implement the change in index at 5500 A rest-frame. 

Our objects are all relatively bright, and so contamination 
from the host galaxy is expected to be negligible (see & how¬ 
ever AGN spectra are also contaminated across the spectral range 
by broad, blended Fe II multiplets. Modelling of these features is 
often accomplished by making use of an empirical template de¬ 
rived from spectral analysis of the Type I AGN I Zwicky 1. We use 
that of |Veron-Cetty et all (2004). There are two free parameters; 
the width of the convolving Gaussian, and the normalisation. We 
fit the broken power-law and Fe II pseudo-continuum to the spectral 
regions in between the strong emission lines - we use: 4000—4300, 
4400-4750, 5050-5800, 5950-6300 and 6800-7500 A. Our fit¬ 
ted continuum is subtracted to leave the Balmer emission line spec¬ 
trum. 

We then fitted the broad permitted emission lines Ha, H/3 and 
H 7 and the narrow, forbidden lines [O III]. We see two categories 
of object in our sample; three show strong, narrow [O III] lines, 
and the other eight do not. This will be discussed more in Paper II. 
Following common practice, we fit multiple Gaussian components 
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Figure 2. Spectral decomposition of the Balmer lines for J2328+1500. This object shows reasonable S/N 11) across the spectral range, whilst some of the 
other objects in our sample show significantly noisier spectra. Although we only use the FWHM and line luminosity of Hct in making the BH mass estimate, 
we fit the other Balmer lines and the [O III] doublet to more strongly constrain the fit, particularly of the iron emission. There is some evidence of a red wing 
to the H/3 profile. The extent of this and possible explanations will be explored in Paper II. 


to emission lines, which offers a reasonable approximation to the 
line shape, and given the quality of our data, is perfectly adequate 
for our purpose. Lines are fitted as follows: 


i. Ha is fitted with two Gaussian components - one broad and 
one intermediate. These are free in velocity shift and normalisation. 
For the three objects that show strong narrow [O III], a third, narrow 
Gaussian component, of wavelength and velocity width tied to the 
strong member of the [O III] doublet is also included. 

ii. H/3 is fitted with two Gaussians, with wavelengths and ve¬ 
locity widths tied to the corresponding components in Ha and with 
the same amplitude ratio. A narrow component is included if strong 
[O III] is observed. 

iii. H 7 is fitted for completeness with 2 components, tied in am¬ 
plitude ratio, velocity width and wavelength to the corresponding 
components in Ha. 

iv. [O III] is a doublet; each member is fitted with a single Gaus¬ 
sian, tied together in velocity width, and with an amplitude ratio of 
2.98 i jStorey & Zeippen|2000[ l. 


We obtain the BH mass estimate for our objects from the Ha 
line, using the method derived by |Greene & Ho| i ]20()5fr . It is com¬ 
mon to use H/3, as much study has been devoted to calibrating mass 
estimates from H/3 with RM samples, due to the greater availability 
of H/3 in optical spectra (Ha being redshifted to IR wavelengths in 
AGN above 2 ~ 0.3), and also because Ha profiles may be blended 
with narrow forbidden lines [S II] and [Nil]. However, given the 
limited S/N of our data, we prefer to use the ~ 3x stronger Ha 
profile. We note that we do not detect the [S II] doublet in any of 
our objects (even those showing [Olll]), which is supporting evi¬ 
dence that contamination of the Ha profile by [N II] is likely to be 
small. Whilst H/3 mass estimates are made using FWHMh /3 as a 
proxy for the velocity dispersion, and the continuum luminosity at 
5100 A (L 51ooa) as a Proxy f or the BLR size,jGreene & Ho (2005 1 
found strong correlations between FWHMh /3 and FWHMHa and 
between L 5100j ^ an d Ln a , and so used this as the basis for deriva¬ 
tion of their relation between Ha profile and BH mass: 


Mbh = ( 2 . 0 ±g’j) X 10 6 

/ r \ 0.55±0.02 / „,, rTT , , s, 2.06±0.06 

f Lac \ / FWHMhc \ M 

y 10 42 ergs -1 ) \10 3 kms -1 / 

( 2 ) 

Following the method of |Greene & Ho|f2005[ i and others, we 
use the full, broad Ha component to obtain BH mass estimates. 
J 12 a discussed the merits of drawing a distinction between the 
Gaussian components that form the broad profile, defining these 
as ‘broad’ and ‘intermediate’ components, but eventually they de¬ 
termined that the two combined yielded the most reliable BH mass 
estimates. 


3.2 Results 

We use a Levenberg-Marquardt minimisation method throughout. 
An example of the Balmer region spectral decomposition is shown 
in Fig. [2] Measurement errors are estimated using a Monte Carlo 
technique, where different iterations of the data are generated us¬ 
ing the ‘mean’ (measured) flux value and the error on that value. 
For each different iteration, optimal values are recalculated, and 
this is repeated 100 times to estimate the error of each fitted value. 
This procedure is not perfect, as we add noise to the already noise- 
degraded spectral data, but serves as a suitable approximation. The 
fitting errors are not indicative of the error on the resulting mass 
estimate though, which are made in accordance with Equation [2] 
We tabulate the results of the emission line analysis in Table[3] 


4 BOLOMETRIC LUMINOSITY 

With the BH mass estimates in hand, we could make predictions as 
to the bolometric luminosity of each AGN by fitting the broadband 
SED model OPTXAGNF to the optical/IR and X-ray data for each 
object. We will take the bolometric luminosity to mean the total 
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Table 3. Ha line measurements for the sample. For objects common to the 
Shen and Liu (2012) sample, our measurements agree to within la. 


ID 

FWHMhq 
( km s _1 ) 

log 

(L Ha /erg s" 1 ) 

^[OIII] 

log 

(Mbh /Mo) 

1 

5600i200 

44.85i0.02 

1.629 

9.42i0.10 

2 

2940i90 

44.57i0.01 

1.567 

8.68i0.09 

3 

4680±60 

44.86i0.01 

2.192 

9.25i0.10 

4 

3100±100 

44.40i0.04 

1.545 

8.63i0.09 

5 

5100±200 

45.21i0.01 

1.535 

9.53i0.10 

6 

3200±100 

44.53i0.01 

1.577 

8.73i0.09 

7 

2550i80 

44.56i0.01 

1.500 

8.55i0.09 

8 

2460i20 

44.87i0.01 

1.562 

8.68i0.09 

9 

3390i80 

44.94i0.02 

1.623 

9.OliO.09 

10 

7810i80 

44.81i0.01 

1.539 

9.68i0.10 

11 

5000±70 

44.32i0.01 

1.609 

9.02i0.09 


intrinsic luminosity of the nuclear source, excluding stars in the 
host galaxy and any associated reradiation. 

4.1 X-ray Spectrum 

We first fit the X-ray data with an absorbed power-law. We include 
attenuation attributable to both the Milky Way (fixed) and the host 
galaxy (free). This will allow us to verify that the following section, 
in which we fit the |Done et al.| ( |2012} OPTXAGNF model, gives rea¬ 
sonable values. We calculate Milky Way Nu values using the Lei¬ 
den/Argentine/Bonn Survey of Galactic HI i jKalberla et al.|2005) . 

In Table |4]the X-ray exposure times and count values are tab¬ 
ulated, and the fitted parameters are shown together with 90% con¬ 
fidence limits. 


4.2 Broadband SED Fitting 

There are a number of properties that affect the intrinsic AGN SED; 
these are described in |Done et al.[ ( j2012| >. The observed SED is also 
affected by optical/UV extinction and soft X-ray absorption due to 
interstellar dust and photoelectric absorption (respectively) in the 
Milky Way. We correct for the extinction (reddening) in the Milky 
Way as discussed in §2.2[ and as in §4. 1 1 include a photoelectric 
absorption component (WABS) to model the soft X-ray absorption 
by Hydrogen, Helium, etc. 

The AGN host galaxy is assumed to have similar intrinsic pro¬ 
cesses reddening/absorbing its emission. We can also model these 
components by redshifting models for X-ray absorption (ZWABS) 
and extinction (ZREDDEN), albeit with no means of constraining 
these other than the shape of the SED. We therefore also produce a 
second model in which these components are added. 

In this redshift range our data may not be sufficient to con¬ 
strain the properties of the SX. Equally, it is possible that either our 
X-ray data or optical data may sample this part of the SED, depend¬ 
ing on the coronal radius. Since we know empirically (e.g. J12a) 
that the SX ought to be taking up a significant fraction (~ 70%) of 
the Comptonised energy, we fix its properties to reasonable default 
values, detailed below. 

In both scenarios, the constrained parameters are as follows. 
A discussion of specific exceptions follows. 

i. BH Mass, Mbh : as previously described in fpl we constrain 
BH mass using the method of |Greene & Ho| ( |2005| l based on Ha. 

ii. Redshift, z: as measured in the spectral decomposition. 


iii. Distance, r c : we calculate the comoving distance to each 
source from the measured redshift assuming a flat cosmology (with 
Ho = 70kms _1 Mpc _1 , fi M = 0.27 and n A = 0.73). 

iv. BH spin, a*: initially we will constrain spin to be zero, in line 
with the work of J12a, but specific instances where a spinning BH 
is possible or implied will be explored in Paper II. 

v. Electron temperature of SX, kT e \ fixed at a typical value of 
0.2 keV. 

vi. Optical depth of SX, r: fixed at a typical value of 10. 

vii. Fraction of Comptonised component in SX, /sx: fixed at a 
typical value of 0.7. 

viii. Hydrogen column density (Milky Way), AH, ga i: calculated 
using the Leiden/Argentine/Bonn Survey of Galactic HI I jKalberla] 
jet al. |2005 > 

The fitted parameters are as follows. 

i. Mass accretion rate, rh = Lboi/LEdd- 

ii. Coronal radius, r cor : the radius at which the AD energy is 
reprocessed by SX and PLT. 

iii. Radial extent of AD, r ou t: the outer radius, in R g of the AD. 
In some objects this cannot be constrained - in these r out is given 
as a limit. 

iv. Power-law slope, E: The power-law index of the coronal PLT. 

v. Intrinsic Hydrogen column density, A^int (ZWABS): soft X- 
ray attenuation intrinsic to the host galaxy. 

vi. Intrinsic reddening, E (B — V ) (ZREDDEN): redshifted extinc¬ 
tion curve to account for reddening intrinsic to the host galaxy. 

We also produce for each object a final model in which the 
SX normalisation is permitted to vary. This allows us to test the 
hypothesis that this component is constrained in some objects. 

To summarise, the models fitted are: 

i. Model 1: SX fixed, no intrinsic attenuation 

ii. Model 2: SX fixed, inch intrinsic attenuation 

iii. Model 3: SX free, inch intrinsic attenuation 

Defining the data to be fit in XSPEC is not straightforward. The 
optical/IR spectral data is, as previously mentioned, contaminated 
by emission features, including complex, blended Fe II emission, 
and the Balmer continuum. We therefore selected, using the Van- 
jden Berk et al.| {200l I quasar template, spectral regions free from 
such emission, and binned these narrow wavelength ranges up into 
well-defined photometry points. The error on each point is defined 
as the flux density standard deviation across the bin. The ranges 
used are (where available): 1300-1350, 1425-1475, 1700-1750, 
2175-2225, 3900-4000, 4150-4250, 5600-5700, 6100-6200, 
6900—7000, unless any of these were unsuitable, e.g. for reasons 
of poor S/N. We then fitted the full energy range of available X-ray 
data from XMM EPIC. 

The bluest bins chosen also do not cover the absorption fea¬ 
tures observed by jKaastra et al.||2014j l, which may be expected in 
AGN with high Ah columns (see Table[4]l. We do not tie the AH.mt 
and intrinsic E (B — V) values together. 

The UV photometry from XMM OM and GALEX was not in¬ 
cluded in the modelling process. At the redshifts we are consider¬ 
ing, the UV filters on these observatories cover a broad wavelength 
range over the strong Lyman-a A1216 A emission line, and the 
Ly-a forest beyond, and are therefore not a good indication of the 
continuum level. With high-resolution UV spectral data, it would 
be possible to interpolate over the narrow absorption features in the 
forest, and recover the underlying continuum (e.g. jFinn et al.|20l4) > 
but this would require observations with e.g. HST/COS. 
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Table 4. The model properties for the X-ray spectrum. We fit absorption components for both the Milky Way (fixed) and the host galaxy (free). In many objects 
the host absorption is poorly constrained, due to the high redshift (only the tail of the absortion profile is sampled) and limited number of counts. The EPIC 
count errors are 1 <x, and the errors on the model parameters are the 90% confidence limits, in line with convention in X-ray astronomy. For similar reasons, 
we also quote I' and its uncertainty to two decimal places. 


ID 

Exp. Time 

XMM EPIC Cts 

Nh, MW 

Nh, int 

r 


(s) 


(xl0 2 ° cm- 3 ) 

(xlO 20 cm" 3 ) 


1 

12519 

410±20 

2.64 

Q + 16 

u -0 

? 34+0-24 
z - jh -0.22 

2 

21 304 

850±30 

1.83 

0+ 11 

-0 

2 56 +0 ' 19 
Z - JO -0.16 

3 

10 523 

142±17 

3.89 

50 +11 ° 

^-50 

2 Q9+0-79 
-0.52 

4 

7 179 

203±18 

2.58 

11+ 18 

il -ll 

2 07 +0,12 
z - u/ -0.13 

5 

15 781 

7200zb 100 

4.48 

45+ 13 

1 QQ+0-09 
a -^-0.08 

6 

19 591 

550±30 

4.04 

20 +40 

ZU -20 

9 99+0.41 

7 

154 071 

4950±80 

1.73 

0+5 

-0 

2-30t°;°o^ 

8 

8 117 

760±30 

1.31 

0+!9 

U -0 


9 

23 543 

1130±40 

1.24 

0 +!7 

U -0 

2 20 + °' 19 
z - zu - 0.12 

10 

124 956 

1350T50 

3.85 

0+ 13 

u -0 

1 44+0- 12 

A -^-o.n 

11 

34 705 

710±30 

4.00 

0+ 19 

u -0 

2 iq+0.12 

z ‘ iV -0.13 


Wed 


1.29 

0.89 

1.71 

0.60 

1.32 

0.77 

1.54 

1.45 

0.84 

1.14 

0.62 


4.3 Results 

The results of the SED fitting procedure are as follows. In Table[5] 
the best-fitting parameters for each model are tabulated, including 
the mass accretion rate in yr -1 . In Table[6]the key properties of 
these SED models are listed. In the manner of J12, we have calcu¬ 
lated K 2 —10 keV — ^bol/^- / 2 —10 keV and ^ 5 iQQj( — Tbol/TAiooA- 

the 2 — 10 keV and 5100 A bolometric correction coefficients, as 
these are commonly used proxies for the bolometric luminosity. We 
also give qox (e.g. Lusso et al.|2010| >. The uncertainties quoted are 
the 90% confidence limits, as is conventional in X-ray astronomy. 
We have estimated these using the Fisher matrix, which gives an 
indication of the measurement error. It should be remembered that 
this does not take account of the systematic errors, to which the 
main contributors will be the uncertainties on the mass estimate 
and the flux calibration. 

We note that some of the models show the SX component to 
be unconstrained by the data, and in five objects this is manifested 
by large fitting errors on r cor and /gx in Model 3. For Model 3, in 
the interest of limiting the number of free parameters, and allowing 
the model to converge to a meaningful minimum, we lock IVh, mt, 
and r 0 ut to the Model 2 values. We make an exception to this rule 
for J0839+5754, as the X-ray data is sufficient to well-constrain 

Afit.int- 

All the SED models are plotted, together with the observa¬ 
tional data, in Appendices[A][B]and[C] 

We see a rather limited range of SED shapes, with all but one 
object being disc-dominated, similar to the SEDs of NLSls in the 
J12 sample. The lowest mass objects, J0043+0114, J0157—0048, 
J1021+1315, J1044+2128, J1240+4740 and J1350+2652 have 
unsampled SED peaks. The red wing of the accretion disc is 
better constrained than for the J12 objects though, owing to the 
lower fraction of host galaxy contribution in these high luminosity 
quasars. 

The objects with BH masses > 1O 9 M 0 - J0041—0947, 
JO 118-0052, J0839+5754, J2328+1500 and J2332+0000 - all 
have observational data extending close to or at their SED peaks, 
enabling reliable estimates of bolometric luminosity in these ob¬ 
jects. 

In many of the objects in our sample, it can be seen that some 
combination of host galaxy contribution and dust reradiation be¬ 
come significant redwards of Ha. The hot toroidal dust component 


will be studied in greater detail in Paper II. Further discussion of 
the host contribution is presented in §5.2| 


5 DISCUSSION 

5.1 Model Suitability and Implications 

We find that the |Done et al. | ( [201 2ft model is able to fit the IR to 
X-ray continuum of our sample of 11 1.5 < 2 < 2.2 AGN. We 
agree with the results of |Capellupo et al.| l |2015| in that many of the 
objects can be modelled in the optical-IR regime by a geometrically 
thin, optically thick AD. In eight objects, constraints are put on the 
outer disc radius, which are compatible with considerations of the 
radius at which self-gravity truncates the AD < |Laor & Netzer| 1989) . 
We note that the presence of a SX, observed and characterised by 
studies of local AGN, is both more physical and necessary to better 
define this continuum. The properties of this SX are related to the 
total energy of the Comptonised component, and modelling the X- 
ray spectrum in addition to optical/IR data is important to infer 
information about the SX. 

Using our SED model we are able to place useful constraints 
on the bolometric luminosity for at least five of the objects in our 
sample. We believe that considerable uncertainties may arise if one 
assumes that the mass accretion rate is adequately estimated simply 
by the use of bolometric correction coefficients jCapellupo et al.| 
|2015| >, as we infer a large spread in those parameters within our 
sample. In Table[5]we show K 2-10 keV and ^siooA’ two commonly- 
used proxies for the bolometric luminosity, for the sample. Though 
our sample is not large, we see a large range of values in all three 
models - around a factor of 10 between the minimum and maxi¬ 
mum. If we only consider the five objects with constrained SED 
peaks, this range is a factor of two in the k 5100j ^, and a factor 
of 10 in ft 2 _iokeV> in spite of the similar masses/accretion rates 
of these five AGN. This echoes the findings of |Elvis et al.| ( fl994j ) 
(who give the similar Lboi/T 2 500 A factor) and J12, and suggests 
that the spread is larger than the ~ 20% stated in Capellupo et al. 
< 12015} . This cannot be solely due to the ~ 0.1 dex error on our 
mass estimate. We therefore suggest that BH spin is not the only 
property that cannot be estimated from singular properties of the 
optical spectra. 

Our Model 3 provides the best fit to the data in all objects, 
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Table 5. The optimum fitted parameters for the various SED models. Uncertainties quoted are the 90% confidence limits, as is conventional in X-ray astronomy, 
and are estimated using the Fisher matrix. As such, they are only indicative of the true measurement error. 


ID 

Ah, int 

(10 22 cm -2 ) 

E (B - V) 
(mag) 

m = 

Lbol /^Edd 

M 

(M 0 yr- 1 

^cor 

) (R g ) 

Pout 

(R s ) 

r 

/sx 

-^reduced 

Model 1: No intrinsic attenuation, SX fixed (SX parameters: kT e = 1 

3.2 keV, r = 10, / S x = 0.7) 




1 

(0.0) 

(0.0) 

0.389±0.018 

44±2 

23±3 

>1000 

2.21±0.03 

(0.7) 

10.4 

2 

(0.0) 

(0.0) 

2.79±0.05 

56.8±1.0 

11 .6±0.6 

790±80 

2.53±0.09 

(0.7) 

0.97 

3 

(0.0) 

(0.0) 

0.39±0.06 

30±5 

24±9 

300±60 

2.45±0.03 

(0.7) 

2.39 

4 

(0.0) 

(0.0) 

2.22 ±0.07 

40.3±1.4 

9.87±0.16 

330±30 

1.94±0.05 

(0.7) 

3.61 

5 

(0.0) 

(0.0) 

0.231±0.018 

33±3 

88.3± 1.4 

> 1000 

1.86±0.04 

(0.7) 

3.50 

6 

(0.0) 

(0.0) 

1.24±0.03 

28.3±0.6 

13.2±0.6 

1800±600 

2.14±0.11 

(0.7) 

2.60 

7 

(0.0) 

(0.0) 

2.9±0.2 

44±3 

10.2±0.4 

>10000 

2.27±0.03 

(0.7) 

2.63 

8 

(0.0) 

(0.0) 

2.10±0.04 

42.7±0.8 

15.9± 1.1 

>10000 

1.79±0.07 

(0.7) 

6.38 

9 

(0.0) 

(0.0) 

1.410±0.016 

60.1±0.7 

10.8±0.3 

700±60 

2.18±0.07 

(0.7) 

1.78 

10 

(0.0) 

(0.0) 

0.0371±0.0011 

7.6±0.2 

21 ±3 

56.0±1.9 

1.63±0.06 

(0.7) 

1.77 

11 

(0.0) 

(0.0) 

0.317±0.010 

14.2±0.5 

20±3 

237±14 

2.18±0.08 

(0.7) 

0.89 

Model 2: Incl. intrinsic attenuation, SX fixed (SX parameters: fcT e = 

= 0.2 keV, r = 

10, /sx = 0.7) 




1 

0.0±0.3 

0.051±0.006 

0.61±0.12 

68±14 

25±9 

240±40 

2.38±0.03 

(0.7) 

5.01 

2 

0.0±0.3 

0.015±0.010 

3.3±0.4 

67±8 

10.9± 1.3 

790±60 

2.52±0.19 

(0.7) 

0.95 

3 

0.15±0.12 

0.025±0.015 

0.50±0.14 

38± 11 

25±11 

260±60 

2.50±0.06 

(0.7) 

2.48 

4 

0 .20±0.10 

0.060±0.015 

4.1 ±0.7 

75± 12 

9.0±0.3 

400±40 

2.07±0.07 

(0.7) 

0.51 

5 

0.18±0.06 

0.065±0.003 

0.338±0.007 

48.8±1.0 

80.1±1.2 

>1000 

1.86±0.05 

(0.7) 

1.45 

6 

0 .2±0.2 

0.024±0.010 

1,58±0.16 

36±4 

13±3 

970±130 

2.32±0.22 

(0.7) 

2.52 

7 

O.OOiO.ll 

0.033±0.006 

4.1 ±0.3 

62±5 

9.4±0.2 

>10000 

2.26±0.06 

(0.7) 

2.54 

8 

0.01±0.17 

0.052±0.014 

3.06±0.12 

62±2 

13±2 

>10000 

1.80±0.13 

(0.7) 

1.77 

9 

0.0±0.3 

0.030±0.008 

1.98±0.18 

84±7 

9.8±0.5 

500±30 

2.18±0.13 

(0.7) 

1.62 

10 

0.0±0.3 

0.094±0.015 

0.067±0.007 

13.6± 1.4 

19±4 

53±2 

1.50±0.09 

(0.7) 

1.30 

11 

0.0±0.3 

0.023±0.015 

0.42±0.08 

19±3 

16±4 

211± 17 

2.18±0.14 

(0.7) 

0.88 

Model 3: Incl. intrinsic attenuation, SX free (SX parameters: kT e = 

0.2 keV, r = 10, /sx — free) 




1 

(0.0) 

0.051 ±0.005 

0.60±0.04 

67±4 

25±2 

(240) 

2.17±0.21 

0.83±0.09 

4.69 

2 

(0.0) 

0.011±0.008 

3.2±0.3 

64±6 

15±40* 

(790) 

2.48±0.14 

0.9±0.3* 

0.93 

3 

(0.15) 

0.028±0.013 

0.51 ±0.05 

38±4 

26±5 

(260) 

2.47±0.09 

0.74±0.10 

2.14 

4 

(0.20) 

0.060±0.014 

4.1 ±0.7 

75± 12 

9.0± 1.1 

(400) 

2.07±0.06 

0.7±0.2 

0.42 

5 

0.51±0.07 

0.051 ±0.006 

0.336±0.009 

48.6±1.3 

82.7 ±1.7 

©1000) 

2.06±0.05 

0.56±0.05 

1.29 

6 

(0.2) 

0.024±0.008 

1.58±0.07 

36.1± 1.5 

10±80* 

(970) 

2.32±0.14 

0.7±3* 

2.38 

7 

(0.0) 

0.033±0.006 

4.1 ±0.3 

62±4 

9.9± 1.8 

©10000) 

2.25±0.05 

0 .8±0.2 

2.53 

8 

(0.01) 

0.056±0.004 

3.16±0.10 

64±2 

10 ±10* 

©10000) 

1.80±0.08 

0.7±0.7* 

1.70 

9 

(0.0) 

0.030±0.005 

1.97 ±0.11 

84±5 

10±14* 

(500) 

2.19±0.09 

0.7±2* 

1.58 

10 

(0.0) 

0.091±0.015 

0.066±0.006 

13.5± 1.3 

17 ±4 

(53) 

1.49±0.06 

0.64±0.16 

1.28 

11 

(0.0) 

0.023±0.009 

0.42±0.05 

19±2 

15±40* 

(211) 

2.18±0.08 

0.7±1.2* 

0.85 


*Large error indicative of unconstrained SX parameter 


judging from the Xred fitting statistic, which takes into account the 
increased number of free parameters in Model 3 versus Models 
1 and 2. In some objects the x?ed value is only marginally lower 
than the Model 2 value, which is indicative of a poorly constrained 
SX that does not benefit from the additional parameter freedom. 
Nonetheless, for the benefit of the objects in which the SX compo¬ 
nent is constrained (five objects), the additional freedom in Model 
3 makes this the model of choice. 

Many of the lower mass objects in our sample are predicted 
to have super-Eddington mass accretion rates, akin to the NLSls 
(e.g. J12a) and ULXs (e.g. Sutton et al. |2013| we observe locally. 
However, we have not yet explored high-spin SED models in our 
study; this will be addressed in Paper II. It is possible that there 
may be some model degeneracy between spin and mass accretion 
rate, and that this contributes to the range of spins predicted by 
(Capellupo et al .|2015| ) in their sample. We will therefore explore 
the effect on accretion rate of having higher spin BHs to assess 


whether the super-Eddington rates we have thus far predicted do 
indeed make these objects high mass NLS1 analogues, or whether 
it is more likely that the additional energy arises from moderately 
or highly spinning BHs. 

Another limitation we have not yet explored is the reliability 
of our mass estimate. It is known that uncertainties on virial BH 
mass estimates are large (~ 0.1 dex or greater) and in our study 
so far we have fixed it at the mean value. Allowing this to vary 
by 1-2 a may well improve the fit, or may add another source of 
degeneracy. Again, this will be explored in Paper II. 


5.2 Host Galaxy Contribution to the Optical/IR Continuum 

Throughout this study we have made the assumption that any con¬ 
tribution to the SED from stars in the AGN host galaxy is likely to 
be negligible. This is a common assumption for typical quasars at 
2 > 0.5 (e.g.|Shen et al.pOll). However, we can test the valid- 
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Table 6. The key properties of the various SED models, including bolometric correction coefficients. 


ID 

l°g(Tbol) l°g(T2-10kev) 

(l°g(erg s 1 )) 

re 2-10keV ^°g(^^2500A) l°g( I '^2kev) 

(l°g(erg s- 1 )) 

“OX 

l°g(AT 510 cbi) 

(logtergs^ 1 )) 

K 5100A 

Model 1: No intrinsic attenuation, SX fixed (SX parameters: / 

fcT e = 0.2 keV, r 

= 10, /sx = 

0.7) 



1 

47.17+0.02 

45.28 

76.7 

46.70 

45.15 

1.60 

46.43 

5.42 

2 

47.271+0.007 

44.86 

257 

46.36 

44.83 

1.59 

45.99 

18.9 

3 

46.98+0.07 

44.78 

159 

46.49 

44.72 

1.68 

46.16 

6.56 

4 

47.103+0.015 

44.78 

212 

46.15 

44.55 

1.62 

45.62 

30.3 

5 

47.07+0.03 

45.69 

24.0 

46.68 

45.43 

1.48 

46.56 

3.20 

6 

46.977+0.009 

44.94 

109 

46.18 

44.78 

1.54 

45.89 

12.1 

7 

47.17+0.03 

44.78 

242 

46.20 

44.67 

1.59 

45.91 

18.0 

8 

47.163+0.008 

45.31 

71.8 

46.27 

45.02 

1.48 

45.99 

15.0 

9 

47.293+0.005 

44.99 

202 

46.57 

44.85 

1.66 

46.26 

10.9 

10 

46.255+0.013 

44.60 

45.3 

45.97 

44.25 

1.66 

45.73 

3.38 

11 

46.646+0.014 

44.79 

72.4 

46.11 

44.65 

1.56 

45.73 

8.16 

Model 2: Incl. intrinsic attenuation, SX fixed (SX parameters 

: fcTe = 0.2 keV, r = 10, / S x 

= 0.7) 



1 

47.33+0.09 

45.27 

116 

46.85 

45.19 

1.63 

46.48 

7.15 

2 

47.36+0.05 

44.86 

312 

46.41 

44.83 

1.61 

46.02 

21.7 

3 

47.08+0.12 

44.82 

181 

46.56 

44.78 

1.69 

46.20 

7.56 

4 

47.38+0.07 

44.86 

332 

46.34 

44.68 

1.64 

45.80 

37.7 

5 

47.233+0.009 

45.84 

24.6 

46.86 

45.58 

1.49 

46.65 

3.80 

6 

47.08+0.05 

44.95 

135 

46.25 

44.85 

1.54 

45.94 

13.7 

7 

47.32+0.03 

44.78 

348 

46.31 

44.66 

1.63 

46.02 

20.1 

8 

47.323+0.017 

45.30 

106 

46.39 

45.01 

1.53 

46.10 

16.7 

9 

47.43+0.04 

44.99 

278 

46.67 

44.85 

1.70 

46.29 

13.8 

10 

46.50+0.05 

44.63 

73.0 

46.22 

44.23 

1.76 

45.86 

4.32 

11 

46.76+0.08 

44.79 

93.5 

46.19 

44.65 

1.59 

45.77 

9.77 

Model 3: Incl. intrinsic attenuation, SX free (SX parameters: 

kT e = 0.2 keV, r 

= 10 , /sx = 

= free) 



1 

47.33+0.03 

45.30 

106 

46.84 

45.16 

1.65 

46.48 

7.12 

2 

47.33+0.04 

44.87 

290 

46.39 

44.82 

1.60 

46.02 

20.2 

3 

47.08+0.04 

44.82 

185 

46.57 

44.77 

1.69 

46.20 

7.62 

4 

47.38+0.07 

44.86 

332 

46.34 

44.68 

1.64 

45.80 

37.7 

5 

47.232+0.011 

45.89 

22.2 

46.81 

45.70 

1.43 

46.63 

4.01 

6 

47.078+0.018 

44.95 

136 

46.25 

44.85 

1.54 

45.94 

13.7 

7 

47.32+0.03 

44.78 

348 

46.31 

44.66 

1.63 

46.02 

20.1 

8 

47.337+0.014 

45.29 

110 

46.40 

45.01 

1.53 

46.11 

16.9 

9 

47.43+0.03 

44.99 

278 

46.67 

44.85 

1.70 

46.29 

13.8 

10 

46.49+0.04 

44.64 

71.4 

46.21 

44.23 

1.76 

45.86 

4.28 

11 

46.76+0.05 

44.79 

93.4 

46.19 

44.65 

1.59 

45.77 

9.76 


ity of this assumption by superposing galaxy SED templates onto 
our faintest source, where the fractional stellar contribution will be 
largest. It is likely that the large galaxies that host the quasars in our 
sample are giant ellipticals, but it is known that starburst galaxies 
have significant energy output in the UV regime, and so we apply 
templates for both of these cases. 

We use two of the galaxy templates of |Polletta et al.| |2007fr 
- that of a 5 Gyr-old elliptical (appropriate for our redshift range) 
and that of the starburst galaxy M82, redshifted as appropriate. In 
terms of normalising these galaxy SEDs, we first assess the great¬ 
est possible contribution in 12328+1500 using the Mbh - ibuige 
relation as presented in e.g. Marconi & Hunt (20031 and DeGraf 
|et al.| ( [2(jl4| >. We test J2328+1500 as it has the highest Mbh, yet is 
our faintest source, and will thus almost certainly show the great¬ 
est contribution to the total SED by the host galaxy. The relations 
of both |Marconi & HuntH2003) and |DeGraf et al.| < |2014| ) predict a 
host galaxy of My — —25. 

We can put an upper limit on the host galaxy contribution us¬ 
ing our SED model and the data. This greatest possible host con¬ 


tribution is shown in Fig. [3] and corresponds to a host galaxy of 
My — —23.3, around 1.7 magnitudes fainter than that predicted 
by the Mbh - Tbuige relation. If this were the case, the contribu¬ 
tion to the host galaxy at the SED peak would be < 2%, even for 
the case of a starburst galaxy (this template is ~ 15 times more 
luminous than M82). An elliptical host would make a negligible 
contribution at the SED peak. 


This result, while representing an extreme case for this ob¬ 
ject, suggests that a host galaxy component may need to be in¬ 
cluded when we model the dusty torus component (evident in the 
WISE photometry) in Paper II. However, the contribution by the 
host galaxy to the total nuclear SED energy is small. For our other 
sources, which are brighter and ought to originate in smaller host 
galaxies (via the Mbh - Tbuige relation), the effect of the host will 
be smaller. 
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Figure 3. Comparison of two host galaxy SED templates against the quasar 
SED for the faintest of our objects, J2328+1500. Host templates have been 
normalised to a V-band absolute magnitude of My ~ —23.3, fainter than 
that predicted from the Mbh — Tbulge relation, but at the maximum pos¬ 
sible contribution permitted by the data (dashed red line). Such a situation 
would imply some red contribution to the total flux of the source, but only 
a small contribution at the SED peak. 


5.3 Variability 

AGN are known to exhibit variability across all wavelength ranges. 
Our study requires that the variability between the optical, IR and 
X-ray observations is not large. A discussion of our approach to de¬ 
tecting and correcting differences between the optical and IR spec¬ 
tral fluxes is given generally in j ]2.2| To summarise, we only see 
evidence for a notable difference between optical and IR flux lev¬ 
els in J0041—0947. The origin of this change may be related to 
the Balmer continuum (as modelled by |Shen & Liu|20T2j l, or poor 
quality photometry of 2MASS (to which the TSPEC spectrum was 
normalised). 

Only a sub-set of properties of the AGN can change over 
timescales of a few years. The BH mass and spin are fixed, and 
changes in the mass accretion rate cannot occur faster than the vis¬ 
cous timescale, which is of the order of thousands of years. Another 
possible source of intrinsic variability in AGN may be tidal disrup¬ 
tion events, in which tidal forces on a star passing close to the BH 
can produce large variations in the observed energy output. 

Extrinsic effects can, in principle, modify the observed SED. 
Gravitational microlensing events by a star in a foreground galaxy 
may affect the total observed flux from an AGN, on timescales of 
months, though such events are predicted to be relatively rare. Cur¬ 
rent models of the toms suggest that it is likely to be clumpy, and 
so a change in the optical depth could occur if a clump were to drift 
into our line of sight. Indeed, major changes in the X-ray column 
density have been observed in several nearby AGN, on timescales 
of months to years (e.g. [Puccetti et a l.| 2007)1Walton et al.|2014fr . 

The properties of the intervening material in the AGN host 
galaxy will also have distinct modifying effects on the observed 
SED. In our SED models, we model attenuation as dust reddening 
of UV/optical/IR spectra, and photoelectric absorption of soft X- 
rays. This has the same wavelength dependence as that of the Milky 


Way, but is redshifted accordingly and has normalisation as a free 
parameter. 

These effects will be explored in Paper II. To display the na¬ 
ture of variability in all our sources, we have plotted the available 
multi-epoch spectral data and photometry in the optical/IR bands 
for each object in Appendix[D] This includes one epoch of photom¬ 
etry from each of the main surveys: SDSS photometric, UKIDSS 
and 2MASS, and all available epochs of spectral data. We supple¬ 
ment this with UV photometry from XMM OM and GALEX all-sky 
imaging survey (AIS), to highlight the uncertainty in these values, 
which as mentioned in ; |2.3| are unreliable due to absorption and 
the presence of emission features. Observation dates for the vari¬ 
ous data sources are tabulated in AppendixjE] We have multi-epoch 
X-ray data for two of our objects - J0839+5754 and J1044—2128. 
Treating each observation as a separate data set, we see no statis¬ 
tically significant evidence for variability in the X-ray spectra of 
these objects. 


5.4 Intrinsic Reddening 


In our analysis, we have made the assumption that host galaxy dust 
extinction (intrinsic reddening) occurs via a similar process to ex¬ 
tinction in the Milky Way. We thus use a redshifted Milky Way 
|Cardelli et al.|p989) extinction curve, which produces apparently 
good reddening correction in all objects, except for J1044+2128. 
We have thus also tested two alternative models for dust extinc¬ 
tion, those of the Large Magellanic Cloud (LMC) and the Small 
Magellanic Cloud (SMC). It is immediately apparent that a better 
continuum fit for J1044+2128 is achieved with the SMC extinc¬ 
tion model, and this is corroborated by the y 2 ed fitting statistic. A 
comparison of different reddening curves is shown for this object in 
Fig.Q There is no evidence for the 2175 A feature in J1044+2128. 

Judging by x? e d only, six objects are best fit with a Milky 
Way extinction curve (an example is given in Fig. [5j, one object 
with an SMC curve, and four objects with an LMC curve. In ob¬ 
jects where the inferred intrinsic reddenin g is small, the difference 
between these y 2 values is marginal. Capellupo et al. (20151 came 
to a similar conclusion that different extinction curves are seen in 
different AGN, although they did not test the LMC model. We thus 
propose that Model 3 can be further augmented by including al¬ 
ternative reddening curves to the Cardelli curve used thus far. We 
will start Paper II by remodelling each SED with the best-fitting 
extinction curve. 

Our only means of constraining the intrinsic reddening is 
the continuum shape, which is certainly a limitation - we see in 
J2328+1500 that an E {B — V) of less than 0.1 mag corresponds to 
an increase in Lboi of 70%. This is a limitation for all such stud¬ 
ies. By investigating the effect of changing our model mass within 
the confidence limits of the mass estimate, as discussed in we 
ought to be able to assess the objects in which changes to the SED 
slope due to reddening are degenerate with small changes in the 
mass estimate. Unfortunately these are most likely to occur in the 
objects with a sampled SED peak. Spin degeneracy may also prove 
to be a contributor. The best approach to test this is studying a larger 
sample in which correlations between, e.g. intrinsic reddening and 
Lboi are directly testable, which could help corroborate or rule out 
such degeneracies and sources of systematic errors. 
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v (Hz) 

Figure 4. Two different extinction curves - the Milky Way and the SMC - 
applied to J1044+2128. The Milky Way curve that we assume throughout 
this work clearly produces an inferior fit to the SMC reddening curve. The 
orange line shows the best-fit SED template, once corrected for intrinsic 
reddening, with the dotted orange line showing the intrinsically reddened 
SED. Similarly, in grey is the reddened optical/IR spectral data, and black 
the dereddened data. This data has been convolved with a 20-pixel Gaussian 
to smooth the data. The Milky Way model prioritises the higher S/N optical 
spectra, as these have a bigger effect on the X^ed. parameter than the 
noisier IR spectra. 



v (Hz) 

Figure 5. An equivalent plot to Fig. ^]for J1350+2652. Here, to highlight 
the difference in continuum shapes implied by the two reddening curves, 
we have fixed both to the same E(B — V) value (0.03 mag). In this object, 
the Milky Way reddening curve produces a noticeably better correction to 
the continuum shape. 


6 SUMMARY AND CONCLUSIONS 

In this study, we have selected a sample of 11 quasars at 1.5 < £ < 
2.2. These objects all have optical, IR and X-ray spectral data, and 
UV photometry. We estimate the BH masses in each object using 
the Ha line profile and method of |Greene & Ho| l j2005} , and then 
fit the energy-conserving, three component SED model one 
|et al.|p012fr to each object. We fit three iterations of this model, 
adjusting the number of free parameters between each one. At this 
redshift range we would expect to observe the peak of the SED, 
due to the both the redshift, and the cooler accretion disc compared 
with AGN with lower mass BHs. 

Our main conclusions are the following: 

i. We observe the SED peak, or close to it, in five objects. We 
find that Model 3, which includes intrinsic attenuation and free SX 
normalisation fits best, allowing for the additional free parameters. 

ii. When used in conjunction with the effects of dust reddening, 
we can accurately model the underlying optical-IR (rest UV- 
optical) continuum, and well-constrain the outer disc radius in 
eight objects. 

iii. In the AGN with lower BH masses, we do not observe the 
SED peak, and in these cases the SX is therefore completely 
unconstrained. As a consequence of this, the model x? e d fitting 
parameters do not differ between the models with SX free and fixed 
for these objects. However, the SX contribution appears to be con¬ 
strained to a varying degree in the five objects with data at the peak. 

iv. Using template SEDs for both luminous elliptical galaxies and 
starbursts we show that the host galaxy contribution is insignificant 
at near to the peak of the SED, but it could contribute a fraction of 
the flux observed redward of Ha. It is very likely that a dusty torus 
also contributes to flux here, judging from WISE photometry. We 
will therefore model this component accordingly in Paper II. 

v. We show that UV photometry alone is insufficient to constrain 
the continuum. Ideally UV spectroscopy, e.g. HST/COS, could be 
used to overcome the uncertainty of Ly-a forest absorption. 

vi. The AGN in our sample generally have high Eddington ratios. 
In this respect they resemble the NLS Is, studied in nearby samples. 
This is expected, as in this redshift range, we preferentially observe 
the brightest AGN, which have high accretion rates. In Paper II we 
will test high-spin SED models, and there may be degeneracy in 
some of our objects between spin and mass accretion rate. 

vii. We identify a range of properties in the best-fitting dust 
reddening component, with SMC/LMC reddening laws providing 
better fits than the Milky Way law in five objects. 

viii. Our analysis provides more reliable estimates of the bolomet- 
ric luminosity, as it uses data from across a large range of wave¬ 
lengths, and utilises an energy-conserving SED model. We highight 
the problems of using a single parameter proxy, such as «2-iokeV, 
as a means to derive Lboi, as we see a large spread in such proxies, 
even in our small sample. We note that the six lowest mass objects 
have unsampled SED peaks, and therefore more poorly constrained 
bolometric luminosities. Having demonstrated the principle of ap¬ 
plying our model successfully to multi-wavelength data, a much 
larger sample will be studied in Paper III to search for relationships 
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between the overall SED characteristics and other specific emission 
line and continuum components. 
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APPENDIX A: SEDS, MODEL 1: SX FIXED, NO 
INTRINSIC ATTENUATION 



z/(Hz) 


Figure Al: Data and SED models for the sample. Here we use Model 1 (SX fixed, without intrinsic attenuation). We also plot the full IR-optical spectrum 
for each object. This spectral data is smoothed for clarity by convolving with a 20-pixel Gaussian. The different SED components are shown using the same 
colour scheme as in Fig.^J The attenuated profile is shown by the dotted grey line. 
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APPENDIX B: SEDS, MODEL 2: SX FIXED, INCL. 
INTRINSIC ATTENUATION 



Figure B1: Data and SED models for the sample. Here we use Model 2 (SX fixed, including intrinsic attenuation). The spectral data is smoothed as in Fig. 
I A1 1 and the same colour scheme is used. 
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APPENDIX C: SEDS, MODEL 3: SX FREE, INCL. 
INTRINSIC ATTENUATION 



Figure Cl: Data and SED models for the sample. Here we use Model 3 (SX normalisation free, including intrinsic attenuation). The spectral data is smoothed 
as in Fig.lAll and the same colour scheme is used. 
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APPENDIX D: SPECTRAL AND PHOTOMETRY PLOTS 
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Figure Dl: All available spectral data plotted with photometry from large surveys. We also show the best fitting (Model 3) attenuated SED profile. At 912 A 
rest frame the photoelectric absorption component cuts the transmitted SED flux to zero. It can be seen in a few objects that have been observed on multiple 
occasions by SDSS/BOSS that variability or inconsistent flux calibration has occurred between observations. Observation dates are given in Appendix [E] 
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Figure D2: Data plots continued. 
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Figure D3: Data plots continued. 
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Table El. Multi-epoch observation dates for all of the objects in our sample. We have searched all large area surveys offering good-quality data. 


ID 


IR 



Optical 



uv 



X-ray 



Instrument 

Type 

Date 

Instrument 

Type 

Date 

Instrument 

Type 

Date 

Instrument 

Type 

Date 

1 

2MASS 

Phot 

1998-10-02 

SDSS 

Phot 

2000-09-25 

GAT.F.X 

Phot 

2006-10-31 

XMM EPIC 

Spec 

2002-01-07 


TSPEC 

Spec 

2010 -01-02 

SDSS 

Spec 

2001-09-10 

GAT.F.X 

Phot 

2006-11-21 





TSPEC 

Spec 

2010-11-28 




GAT FX 

Phot 

2011-03-18 




2 

2MASS 

Phot 

2000-11-29 

SDSS 

Spec 

2000-09-07 

GAT FX 

Phot 

2003-09-16 

XMM EPIC 

Spec 

2010 -01-10 


UKIDSS 

Phot 

2006-11-19 

SDSS 

Phot 

2008-10-02 

GAT FX 

Phot 

2003-09-30 





UKIDSS 

Phot 

2008-11-28 

BOSS 

Spec 

2009-12-21 

GAT FX 

Phot 

2007-10-27 





GNIRS 

Spec 

2013-08-16 

BOSS 

Spec 

2010-09-05 







3 

2 MASS 

Phot 

1998-09-18 

SDSS 

Spec 

2000-09-02 

GAT FX 

Phot 

2008-10-20 

XMM EPIC 

Spec 

2003-07-11 


GNIRS 

Spec 

2004-11-29 

SDSS 

Phot 

2004-09-23 

GAT FX 

Phot 

2008-10-31 





UKIDSS 

Phot 

2006-07-10 




GAT FX 

Phot 

2008-11-16 











GAT FX 

Phot 

2011-10-28 




4 

2MASS 

Phot 

1998-09-29 

SDSS 

Spec 

2000-11-23 

GAT FX 

Phot 

2004-10-11 

XMM EPIC 

Spec 

2005-07-14 


UKIDSS 

Phot 

2005-09-07 

SDSS 

Spec 

2001-09-27 

XMMOM 

Phot 

2005-07-14 





TSPEC 

Spec 

2009-11-07 

SDSS 

Phot 

2003-11-19 

GAT FX 

Phot 

2008-10-20 





TSPEC 

Spec 

2010-11-28 

BOSS 

Spec 

2010-09-10 







5 

2MASS 

Phot 

2000-01-05 

SDSS 

Phot 

2003-10-23 

XMMOM 

Phot 

2006-10-03 

XMM EPIC 

Spec 

2006-10-03 


GNIRS 

Spec 

2013-10-27 

SDSS 

Spec 

2007-11-21 

GAT FX 

Phot 

2007-01-03 

XMM EPIC 

Spec 

2007-04-06 








XMMOM 

Phot 

2007-05-09 

XMM EPIC 

Spec 

2007-05-09 








GAT FX 

Phot 

2010-01-14 




6 

UKIDSS 

Phot 

2010-02-08 

SDSS 

Phot 

2003-01-27 

XMMOM 

Phot 

2003-05-05 

XMM EPIC 

Spec 

2003-05-05 


GNIRS 

Spec 

2014-03-21 

SDSS 

Spec 

2004-02-27 

GAT FX 

Phot 

2006-03-27 











GAT FX 

Phot 

2010-03-14 




7 

2MASS 

Phot 

1998-01-29 

SDSS 

Phot 

2005-03-09 

GAT FX 

Phot 

2006-09-14 

XMM EPIC 

Spec 

2003-05-05 


GNIRS 

Spec 

2014-03-20 

SDSS 

Spec 

2006-12-28 




XMM EPIC 

Spec 

2003-05-28 





BOSS 

Spec 

2012-04-22 




XMM EPIC 

Spec 

2003-12-12 

8 

2MASS 

Phot 

1998-05-16 

SDSS 

Phot 

2003-03-10 

XMMOM 

Phot 

2002 -11-12 

XMM EPIC 

Spec 

2002 -11-12 


TSPEC 

Spec 

2011 -02-22 

SDSS 

Spec 

2004-03-25 

GAT FX 

Phot 

2007-03-04 




9 

2MASS 

Phot 

2000-04-11 

SDSS 

Phot 

2004-06-11 

XMMOM 

Phot 

2004-01-25 

XMM EPIC 

Spec 

2004-01-25 


UKIDSS 

Phot 

2010-03-01 

SDSS 

Spec 

2006-04-23 

GAT FX 

Phot 

2006-04-30 





TSPEC 

Spec 

2011 -02-22 

BOSS 

Spec 

2012-06-27 

GAT FX 

Phot 

2009-05-27 











GAT FX 

Phot 

2011-05-27 




10 

UKIDSS 

Phot 

2007-09-28 

SDSS 

Phot 

2000-09-26 

GAT FX 

Phot 

2004-09-15 

XMM EPIC 

Spec 

2007-12-01 


GNIRS 

Spec 

2013-08-18 

SDSS 

Spec 

2001-11-25 

GAT FX 

Phot 

2006-02-26 











GAT FX 

Phot 

2007-03-28 











GAT FX 

Phot 

2009-09-09 











GAT FX 

Phot 

2009-10-07 




11 

UKIDSS 

Phot 

2006-05-06 

SDSS 

Spec 

2000-10-04 

GAT FX 

Phot 

2004-03-14 

XMM EPIC 

Spec 

2007-12-01 


GNIRS 

Spec 

2013-08-19 

SDSS 

Spec 

2001-10-17 

GAT FX 

Phot 

2006-10-01 








SDSS 

Spec 

2002-09-08 

GAT FX 

Phot 

Many obs 








SDSS 

Phot 

2003-11-19 



in DIS 
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